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We demonstrate genuine three-mode nonlocality based on phase space formalism. A Svetlichny- 
type Bell inequality is formulated in terms of the s-parameterized quasiprobability function. We test 
such tool using exemplary forms of three-mode entangled states, identifying the ideal measurement 
settings required for each state. We thus verify the presence of genuine three-mode nonlocality 
that cannot be reproduced by local or nonlocal hidden variable models between any two out of 
three modes. In our results, GHZ- and W-type nonlocality can be fully discriminated. We also 
study the behavior of genuine tripartite nonlocality under the effects of detection inefficiency and 
dissipation induced by local thermal environments. Our formalism can be useful to test the sharing of 
genuine multipartite quantum correlations among the elements of some interesting physical settings, 
including arrays of trapped ions and intracavity ultracold atoms. 

PACS numbers: 03.65.Ud, 03.65.Ta, 03.67.-a, 42.50.-p 



> 

ov 
o 
in 



X 



Quantum nonlocality is one of the most fundamental 
features of quantum mechanics. It refers to the correla- 
tions that cannot be explained by local hidden-variable 
models that satisfy a set of constraints epitomised by 
so-called Bell inequalities [l| . The violation of a Bell in- 
equality reveals the existence of nonlocality in a given 
quantum mechanical state . 

While originally formulated for bipartite systems, Bell- 
like inequalities have been extended to the multipartite 
scenario, a noticeable example being embodied by the 
well-known inequality proposed by Mermin and Klyshko 
(MK) @. However, the violation of a MK-type inequal- 
ity by a multipartite state does not necessarily imply the 
existence of genuine multipartite nonlocality, as this test 
can be flasificd by nonlocal correlations in any reduction 
of the system's components. In order to demonstrate 
genuine tripartite nonlocality, another type of Bell in- 
equalities should be thus considered such as the one for- 
mulated by Svetlichny Q , which rules out both local and 
nonlocal hidden variable model possibly imposed on any 
sub-parties. It was also noted that a stronger violation 
of the MK type can demonstrate genuine nonlocality for 
the cases with even number of parties Q. Experimen- 
tal demonstrations of genuine multipartite nonlocality 
were firstly achieved by strong violations of a MK in- 
quality with four photon polarization entanglements Q. 
A violation of Svetlichny type inequality was experimen- 
tally demonstrated recently with GHZ-type photon po- 
larization entangled states @. A generalised version of 
Svetlichny inequality was recently proposed and studied 

In this paper, motivated by the growing experimental 
capabilities of controlling and manipulating the state of 
tripartite quantum systems, in the optical lab (ll| and 
beyond, we address the formulation of genuine three- 
mode continuous variable (CV) nonlocality tests in phase 



space. While the phase space provides a natural arena 
for the description of the state of multi-mode CV sys- 
tems, it also provides us with powerful tools for the anal- 
ysis of the quantum correlation features within a given 
state and the possibility to investigate the quantum-to- 
classical transition in a transparent, intuitive way. While, 
historically, the completeness of quantum mechanics has 
been questioned by Einstein, Podolsky and Rosen using 
an argument related to position and momentum of a CV 
system [l2| , the first formulation of a Bell-like nonlocality 
test in phase space was provided by Banaszck and Wod- 
kiewicz using the Wigner and Husimi Q-function (l3| . 
This approach has been generalized and extended in sev- 
eral directions, including the provision of phase-space 
tests valid for high-dimensional systems [14rHa | and the 
addressin g o f general s-parameterized quasiprobability 
functions (l7H19| . Although a considerable body of work 
has been produced to address issues of non locality in 
multi-mode CV states [20l - [26j . a self-consistent formu- 
lation of genuine multimode nonlocality tests in phase 
space is still missing. 

In this paper, we fill this important gap by formalis- 
ing MK and Svetlichny-type Bell inequalities in terms of 
the s-parameterized quasiprobability functions [nj, and 
thus opening the way to the experimental test of multi- 
partite nonlocality of CV systems in their phase space. 
In order to assert the usefulness of our formal tools, we 
study in depth the tripartite nonlocal character of three 
paradigmatic states: a three-mode single-photon entan- 
gled state, a tripartite squeezed vacuum state, and the 
entangled coherent state of three bosonic modes. By ex- 
ploiting the generality of our formulation and exploring 
the properties of these classes of states, we build up com- 
plete tripartite-nonlocality phase diagrams, identifying 
regimes where the presence of genuine three-mode non- 
locality is guaranteed and investigating its character in 
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terms of sharing structure of quantum correlations. We 
complete our analysis by analysing the effects of both lo- 
cal damping and detection inefficiencies on the genuine 
tripartite nonlocal nature of a given state, proving that 
our inequalities are flexible enough to incorporate such 
detrimental effects in a natural way. Although we do not 
explicitly address it here, the extension of our formalism 
to the larger registers of modes is fully straightforward. 
Moreover, as we show in this paper, this investigation 
holds the promises to provide a useful tool for the char- 
acterization of the quantum correlation-sharing structure 
in physical settings of current enormous technological in- 
terests. 

The remainder of this paper is organized as follows. 
Bell-like inequalities for multipartite scenarios are briefly 
reviewed in Sec. |TJ This is then followed by our formula- 
tion of phase-space tests based on generalized quasiprob- 
ability functions, as presented in Sec. |TTJ In Sec. IIIII we 
pass to investigate the violation of such inequality by the 
abovementioned classes of three-mode entangled states, 
whose tripartite nonlocal character is studied against rel- 
evant sources of detrimental effects in Sec. IIVI Sec. Mis 
then devoted to the identification and brief discussion of 
a series of physical setups where our study will be invalu- 
able for the identification of the proper sharing structure 
of quantum correlations. Finally, in Sect. lVll we draw our 
conclusions and highlight potential lines of developments. 

I. GENUINE MULTIPARTITE NONLOCALITY 

In this Section, we discuss the formal approach to the 
derivation of Bell inequalities for the test of multipartite 
nonlocality. In order to fix the ideas, we consider ex- 
plicitly the case of a tripartite system, although the dis- 
cussion here is straightforwardly generalised to the more 
than three subsystems. 

It is instructive to briefly revise the two-party case, 
first. Let us thus assume that two particles, whose 
state can be repeatedly and identically prepared, arc dis- 
tributed to two remote parties, Alice, Bob. Each of them 
performs a measurement chosen out of two: A(A') for Al- 
ice and B(B') for Bob. The corresponding possible out- 
comes are assigned as j (J') with j = a, b, respectively, 
and j G {—1,1}. Under the assumption of the existence 
of local hidden variables, the joint outcome probability 
after many rounds of measurements is written as 

P(ab\AB)= f d\ p(\)P x (a\A)P x (b\B), (1) 

where A is a shared local variable and P\{j\J) is the 
probability that measurement J — A, B has outcome j = 
a, b at a set value of A. The Bell-CHSH parameter can 
then be cast in the form of 

B = a(b + b') + a'(b-b'). (2) 

It is straightforward to check that \B\ < 2. The violation 
of such inequality witnesses the failure of local hidden 



variable theories to describe the correlations between the 
outcomes of the measurements performed by Alice and 
Bob and, thus, the untenability of the assumptions of 
locality and realism. It is well known that if the two dis- 
tributed particles are prepared in an entangled quantum 
state, the Bell-CHSH inequality can be violated by up to 
a factor \f2. 

This reasoning can be extended to the inclusion of a 
third party, say Charlie, who performs measurements C 
and C with outcomes c, d € { — 1, 1} on a third particle. 
In this case, it is possible to combine the outcomes of 
each measurements so as to build the parameter [B[ 

M = \b(c + c>) + \b'(c-c>) (3) 
= a'bc + ab'c + abc' — a'b c , 

which is bounded as \M.\ < 2 according to local hidden 
variable theories. Eq. ([3]) embodies the Mermin-Klyshko 
(MK) inequality, which can be violated by quantum me- 
chanics, thus showing tripartite quantum nonlocality. 

However, any violation of the MK inequality in Eq. J3]) 
does not guarantee the existence of genuine tripartite 
nonlocality: nonlocal correlations between any two par- 
ties out of three are sufficient to exceed the boundary im- 
posed by local hidden variable theories. To avoid this, the 
joint probability for genuine tripartite nonlocality should 
not be reproduced by any local hidden variable assigned 
on the joint measurement of any two out of three parties. 
That is, we define the joint probability 

P(abc\ABC) = f d\J2 V iP(ab)cWP\(ab\AB)P x (c\C)}, 

(4) 

where V performes the cyclic permutation of the triplets 
of indeces (abc) and (ABC), P\(ij\IJ) stands for the 
joint probability of outcomes i,j = a,b,c for the mea- 
surements I, J = A, B, C, and P( a b)c denotes the cor- 
responding probability distribution density. Building on 
this, Svetlichny proposed a combination of measurements 
outcomes of the form Q 

S = Bc±B'c = M±M' 

= a'bc + ab'c + abc' - a'b'c' (5) 
± ab'c' ± a'bc' ± a'b'c T abc. 

Again, it can be proven easily that |5| < 4 under the 
assumptions behind the Eq. (01 . The violation of the 
Svetlichny inequality (JS|) thus signals genuine tripartite 
nonlocality. Quantum mechanics is know to violate the 
bounds imposed by local hidden variables: when the 
probabilities entering Eq. ([5]) are calculated performing 
measurements over a tripartite GHZ or a W state, the 
Svetlichny parameter (S) is larger than 4 (achieving the 
value 4V2 over a GHZ state) . 
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II. PHASE SPACE FORMULATION VIA 
GENERALIZED QUASIPROB ABILITY 
DISTRIBUTIONS 

In this Section we re-formulate multipartite Bell in- 
equalities in terms of generalized quasiprobability func- 
tions defined in the phase space of CV systems. Our 
approach will be based on the formalism introduced in 
Refs. [13, EH and provides the complement to the seminal 
analysis by Banaszek and Wodkiewicz |l3| of bipartite 
CV states. The s-parameterized quasiprobability func- 
tion for a given single-mode state p reads p7j - [l9l | 



W(a;s) 



7r(l - s) 



Tr[pTL(a;s)}, 



(6) 



where fl(a; s) = £)£L [(s + 1 )/( s - l)H a , n){a, n\, and 
\a,n) = D{a)\n) is the n-photon state of a boson de- 
scribed by annihilation and creation operator a and a', 
displaced in phase space by the Weyl operator D(a) = 
expfaat — a*a\. W(a;s) reduces to the P, Wigner, and 
Husimi Q function for s = 1,0,-1, respectively 27 1. 

Suppose that Alice, Bob and Charlie, independently 
choose one of two observables, denoted by A a , Bi,, C c 
respectively, where a, 6, c = 1,2, where no restriction is 
placed on the number of possible outcomes. The local 
measurement operators are written as 





a a 


for (J, j) 


= (A,a) 


Jj = 0{5j] s) with Sj = < 


Pb 


for (J,j) 


= (B,b) 




J- 


for (J,j) 


= (C,c) 



any local observable is bounded by 1 for any s < 0, the 
expectation value of the MK operator given in Eq. (fit))) 
is such that \M\ < 2 for any local hidden variable theory. 
From Eq. ((6]), the expectation value of II(a; s) is propor- 
tional to the s-parameterized quasiprobability function 
as (fl(a; s)) = 7r(l — s)W(a; s)/2. Note that we do not 
consider the case s > when the eigenvalues of n(a; s) 
are unbound. These considerations allow us to calculate 
the explicit form of the s-parameterized MK parameter 
for the two ranges of values of s considered so far. Such 
expressions are too long and uninformative to be given 
here and are deferred to Appendix 1. 

Starting from Eq. ([5]) and the expression in Eq (|A-1[) 
we can easily construct the s-dependent Svetlichny pa- 
rameter such that 



\S S \ = \M S ±M' S \<4. 



(11) 



and the Hcrmitian operator 



This is a generalization of the Svctlichny-likc inequal- 
ity formulated for s-paramcterized quasiprobability func- 
tions. For example, it reduces to the Svetlichny inequal- 
ity formulated for the Wigner function [2(| when s = 
and for the Husimi Q-function for s = — 1. The full form 
of such Svetlichny parameters is given in Appendix 1 . Us- 
ing the MK and Svetlichny inequalities in Eq. (|A-1|) and 
Eq. (fTTj) , we can test genuine three mode nonlocality for 
arbitrary systems represented in phase space. 

In the remainder of the paper we apply our phase-space 
formalism to the study of a few paradigmatic examples 
that will serve as useful benchmarks. 



III. TESTING GENUINE NONLOCALITY FOR 
THREE-MODE ENTANGLED STATES 



I 211(a; s) — 1 it s < — 1 

with 1 the identity operator. The possible measurement 
outcomes of 0(a; s) are given by the eigenvalues 



(l-s)(f^)" + s if-l<s<0, 
2(f±i)"-l if s < — 1, 



and their eigenvectors are the displaced number states. 
The maximum and minimum measurement outcomes of 
0(a; s) for any non-positive s are A max = — A min 
1 . For s = we have 6(a; 0) = n(a;0) = 
Sn°=o( — l)™! 01 ' n )( a > n \i which is the displaced parity 
operator. On the other hand, for s = — 1 we have 
0(a;~- 1) = 2|a)(a| — 1, i.e. a projector onto coherent 
states. The MK and Svetlichny parameters can be con- 
structed using the measurement operators Jj. Indeed, 
from Eq. ([3]), the MK parameter is now 



M = c 112 +c 1 2i+e 211 -c 2 



(10) 



where C a t c = {A a (g> Bf, (g) C c ) is the correlation function 
for measurement outcomes. As the expectation value of 



Here, we investigate genuine tripartite nonlocality us- 
ing our phase-space formulation of Svetlichny-like in- 
equalities by addressing three exemplary cases of three- 
mode entangled photon states. Specifically, we consider 
the three-mode single-photon entangled state 



(9) \4>) = \i- 



7 (|o)|o)|i) + |o)|i)|Q)) 

V2 




|1)|0)|0), (12) 



with p £ [0, 1], which becomes a W-state for p = 1 and a 
two- mode single-photon entangled state for p — 0. The 
state can be generated when a single photon enters a 
three- mode beam splitter (fritter). Its s-parameterized 
quasiprobability function can be written as 



W^(a,/3,7;s) 




§(l-f)0 + 7) 

(13) 

Quantum nonlocality of a two-mode single-photon state 
has been demonstrated theoretically [l?} , but so far there 
have been no addressing of the genuinely nonlocal nature 
of its three-mode counterpart. 
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The second example that we consider is the three-mode 
squeezed vacuum (3MSV) states introduced in Ref. [21 



which can be generated by combining three single-mode 



squeezed states (with identical degree of squeezing r) at 
a tritter [U [29[. Its s-parameterized quasiprobability 
function is 



W 3M sv(ai,P,T, s ) 



7r 3 (l + s 2 - 2s cosh 2r) 3 / 2 



cxp 



3(1 



+ {aj + 01 + 7l 2 - 4a,/3 4 - 4/3 l7l - 4 7l a, 



2s cosh 2r) 



3(s — cosh2r)(|a| 



+ |7| 2 ) 



Pr - Ir + 4av/3 r + 4/3 r7r + 4 7r a r ) sinh 2r 



(14) 



where S r and 6i denote the real and imaginary parts of 
5 = a, /3, 7 . Such Gaussian state has been introduced as 
the CV version of the GHZ entangled state [2l| , although 
for small squeezing it contains some residual two-mode 
entanglement and it can thus be also regarded as a rep- 
resentative of W-class entanglement 22] . 

Finally, we address entangled coherent states (ECS) of 
the form 

|ecs) =^(|c)|C)IC> -I -01 -01-0), (is) 

where £ is the amplitude of a single coherent state (we 
assume for simplicity that ( G E) and Af is a normaliza- 
tion constant. A scheme for the generation of these states 
was introduced in Ref. [28[. Although states |ECS) have 
been regarded as the entangled coherent state version 

it has been also known 
TSJ reduces to Eq. (13]) 



28] 



of GHZ-class entanglement 
that, in the limit of ( — s- 0, Eq. 
with p = 1. Therefore, depending on the amplitude of 
the state components, it might be possible to identify 
two distinct behaviors of the Svetlichny parameter, as- 
sociated with either GHZ or W-class entanglement. The 
s-parametcrized quasiprobability function of |ECS) is 



We CS (a,/3, 7 ;s)= 



8J\f 2 



\l-sf 



-2e z 



\g=± 
4C(cKi + /3i + 7 i) 



i+gcr+li+eO 



1 - 



(16) 

It would also be tempting to study the ECS states 
A/"'(|0) |0) |C) + |0)|O]0> + |C)|0)|0)), sometimes referred 
to as W-type ECS [2g]. However, the non locality of this 
class of states cannot be revealed by relying on Bell-like 
tests based on local displacement in phase space, as con- 
sidered here, and require the use of local operations based 
on Kerr-like nonlinearities instead (28j . 

The three paradigmatic classes of states introduced 
above will be now studied against their tripartite non- 
local character. Later on in this Section we consider 
the behaviour of both MK and Svetlichny-like param- 
eters against the influences of dissipation affecting such 
states, which will allow us to comment on the robustness 
of genuine tripartite nonlocality. 



A. Demonstrating genuine three-mode nonlocality 
in phase space 

We start investigating the violations of both MK and 
Svetlichny type inequalities for each of the entangled 
states discussed above. As in our formalism quantum 
nonlocality is independent of the local measurement set- 
ting associated with parameter s, the violation of a 
Svetlichny-like inequality for a suitably chosen s guaran- 
tees the presence of genuine three-mode nonlocality in the 
specific entangled state being studied. In the following, 
we refer to the genuine tripartite nonlocality exhibited 
by states belonging to the GHZ- and W-class entangled 




FIG. 1: Behavior of the s-parameterized MK and Svetlichny 
parameter \M\ and |<S| for (a) the single-photon W state \(j>) 
(with variable p), (b) the three- mode squeezed vacuum state 
(r being the degree of squeezing), and (c) the GHZ- type ECS 
(£ is the amplitude of each coherent-state component). 
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states as GHZ-type and W-typc nonlocality , respectively. 

The genuine tripartite W-type nonlocal nature of a 
single-photon entangled state |</>) for proper choices of p 
can be appreciated from the right panel of Fig.[IJa): Both 
the MK and Svetlichny inequalities arc violated maxi- 
mally by a test designed using the Husimi Q function 
(i.e. for s = —1). However, while the MK inequality is 
violated for all values of s (Fig. Q] (a), left panel), the 
Svetlichny one is violated only close to s = — 1, implying 
the relevance of the choice of suitable local measurements 
to demonstrate genuine three-mode nonlocality. Clearly, 
when p = the state is solely endowed with two-mode 
entanglement, which yields the violation of the MK test 
only. The Svetlichny parameter depends monotonically 
on p and achieves its maximum value at p = 1 . 

For the three-mode squeezed vacuum states, the 
Svetlichny inequality can be violated only in a region 
of values very close to s = [cf. Fig. QJb), right panel], 
where the dependence of S on the degree of squeezing is 
monotonic and where observables embodied by local dis- 
placed parity operations (as it is the case for phase-space 
formulation of Bell-like tests based on the Wigner func- 
tion) arc optimal. From the observation in Ref. [22j that 
the entanglement of any two-mode subsystem disappears 
for sufficiently large degree of squeezing r, we claim for 
GHZ-type nonlocality, in this case. For small values of 
r, though, W-class entanglement becomes dominant in 
the state. For example, it is known that a three-mode 
squeezed vacuum state with r ~ 0.5 contains W-type en- 
tanglement dominantly as its residual two-mode entan- 
glement becomes maximum (22j. This is in agreement 
with our analysis of the MK parameter M. , whose formu- 
lation in terms of the Q function (s = —1) maximizes the 
degree of violation of the MK inequality, as shown in the 
left panel of Fig. [IJb). However, note that W-type non- 
locality does not appear in this range of r, implying that 
not all W-class entanglement can yield genuine W-type 
nonlocality in phase space. 

For the ECS in Eq. (T5]), both GHZ and W-type non- 
locality can emerge with an appropriate choice of s, as 
shown in Fig. [ljc). For small amplitudes genuine tri- 
partite nonlocality is observed for a test built using the 
Q function. This part of the tripartite non locality phase 
diagram for an ECS is originated by the W-class entan- 
glement of a single-photon state and thus reveals W-type 
nonlocality. Needless to say, the results associated with 
( — > are the same as those obtained for p = 1 in state 
|</>) and shown in Fig. [TJa) (right panel). On the other 
hand, GHZ-type nonlocality is achieved in a very nar- 
row region close to s = 0. In this case, a test based on 
the Wigner function (i.e. the use of local parity mea- 
surements) is optimal for demonstrating CV GHZ-type 
nonlocality. As shown in Fig. [2J the W- and GHZ-type 
nonlocal characters are interchanged at £ ~ 0.455; past 
this point, the W type (best revealed by taking s = — 1, 
i.e. using a test based on the Husimi Q function) dis- 
appears to leave room for a GHZ-type character (to be 
tested by the Wigner function). Therefore, genuine tri- 



partite nonlocality can be observed for all ranges of C 
with suitably chosen local measurement setting. 



IV. EFFECTS OF DETECTION INEFFICIENCY 
AND LOCAL DAMPING ON GENUINE 
NONLOCALITY 

Here we investigate the effect of the two important 
detrimental mechanisms on the genuine three-mode non- 
locality demonstrated in Sec. IIII1 Realistic implementa- 
tions of CV nonlocality tests should take into account 
both the inefficiency of the detectors used to collect the 
statistics and the potential interaction of each local mode 
participating to the state being probed with environmen- 
tal baths. Both effects can be effectively described as 
changes of the parameter s used in our formalism. In 
particular, as it is shown in Appendix 2, the single-mode 
s-paramctcrized quasiprobability function measured by a 
detector with efficiency 7/ £ [0, 1] is given by a 1/rj rcscal- 
ing of the ideal quasiprobability function that should then 
be parameterized by [3(| 



Therefore, the MK- and Svetlichny-type inequalities un- 
der the effect of detection noise can be formulated by 
replacing the quasiprobability distribution in Eq. (|A-1[) 
and Eq. (fTTj) with the measured quasiprobability func- 
tions, ej^ W 3 (a, (3,j;s) -> W 3 (a, /3, 7; s') / (jq a -q b ^ c ) (j = 
1,2,3) with r\ a .b,c the detection efficiency of the detec- 
tors at Alice's, Bob's, and Charlie's site, respectively (the 
marginals being calculated accordingly). For simplicity, 
in what follows we will restrict the attention to the case 
of Va.b.c — V- Similar transformations are enough to take 
into account the effects of the interactions of each lo- 
cal mode with its own thermal reservoir (at equilibrium 

6.0 




12 3 4 



FIG. 2: Genuine tripartite nonlocality phase diagram for 
a three-mode ECS. We plot the Svetlichny parameter |<S| 
against the amplitude of the coherent-state components 
The W-type and GHZ-type nonlocal characters, best revealed 
by tests built on the Husimi (s = — 1, dashed line) and Wigner 
function (s = 0, solid line), respectively, cross for £ ~ 0.455. 
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FIG. 3: Tripartite non locality phase diagram for a three- 
mode ECS against the detection efficiency 77. As in Fig. [2] 
dashed (solid) traits are for s = — 1 (s = 0) 



with h thermal photons) through non-zero temperature 
amplitude-damping mechanisms 31] occurring at rate T 
(assumed to be the same for all modes, for the sake of 
argument). In this scale, we should operate the replace- 
ment 



Ws(a,/3, T ,s) 



1 w ( _a P 2_ '/ \ 

t"(r) W3 \t(rrt(ryt(ry S{T) 



(18) 

with i(r) = \J\ — r 2 (r) = %/ e~~ Tr , r the mode-bath in- 
teraction time and 



s'(r) 



s-r 2 (r)(l+2n) 



Equipped with this tools, we can now address the effects 
of both detection inefficiency and local damping on the 
W- and GHZ-type nonlocality. 

In Fig. [3j we plot the effect induced by imperfect de- 
tection on the behavior of the Svetlichny parameter for a 
three-mode ECS. Linking to our previous considerations, 
we focus on the trade-off between W- and GHZ-type non- 
locality. We use the trends shown in Fig. [2] as a bench- 
mark for the functional behaviour associated with 77 < f. 
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FIG. 4: Behavior of the MK parameter [panel (a)] and the 
Svetlichny one[panel (b)] against the section efficiency r\ for 
a three-mode ECS (£ = 1 .0, O.f ) and a three-mode squeezed 
state (r = 0.5). 



FIG. 5: Dynamics of genuine GHZ-type (£ = 1,2) and W- 
type (C = 0.1) nonlocality of a three-mode ECS interacting 
with individual thermal baths at equilibrium (mean occupa- 
tion number n — 0). 



Clearly, GHZ-type nonlocality is somehow more prone to 
the effects of imperfect detection, as the solid-line traits 
in Fig. [3] (associated with Wigner function-based non- 
locality tests, i.e. for s = 0) disappear quickly for even 
very weakly inefficient detectors. The W-type nonlocal 
character appears to be quite more robust: r) > 0.97 
is required to demonstrate GHZ-type nonlocality for the 
entangled coherent state with ( = 1.0, while W-type non- 
locality for the state with ( = 0.1 can be observed with 
77 > 0.955. Needless to say, this fragility is partly due 
to the large amplitude of the coherent-state components 
needed to test GHZ-type nonlocality. Continuing our 
analysis, we see that, expectedly, in comparison with the 
MK parameter, the Svetlichny one is quite more sensi- 
tive. In Fig. HI we find the threshold in the detection 
efficiency that is needed in order to observe the violation 
of both MK and Svetlichny inequalities by a three-mode 
ECS. The verification of genuine W-type nonlocality with 
a W-class ECS state (( = 0.1) requires 77 > 0.955, while 
T] > 0.78 is enough to observe the violation of the MK 
inequality. These considerations hold, qualitatively, also 
when the effector local damping is considered [cf. Fig. [5] 



Our final confederation in this respect is that the 
robustness shown by the MK parameter for a given 
three-mode state does not always imply an equally ro- 
bust Svetlichny parameter. For example, the three-mode 
squeezed vacuum state with degree of squeezing r = 0.5 
exhibits maximal W-class entanglement [22| and has an 
associated MK parameter that is more resilient to detec- 
tion inefficiencies than the ECS with amplitude C = 0.1 
as shown in Fig. [H However, it does not yield any gen- 
uine W-type nonlocality, even with perfect detectors, in 
contrast with the some how more fragile ECS. 



V. PHYSICAL SETTINGS 

The progresses made in the last ten years in the design, 
manipulation and control of mesoscopic quantum sys- 
tems consisting of hybrid components has now made pos- 
sible the experimental realisation of interfaces between 
devises as diverse as ultra cold atoms and high-finesse 
cavities [32[ or micro mechanical oscillators [33]. On 
the other hand, the development of well-acquired tech- 
niques for the spacial confinement of charge atoms has 
now reached outstanding levels allowing for the quantum 
processing of tens of individual particles [34j | . This is very 
interesting for the purposes of our study. 

Indeed, although throughout this paper we have used 
the language typical of all-optical implementation of 
CV systems, the formal apparatus built in our study 
applies, needless to say, to any effective bosonic sys- 
tem, regardless of its physical embodiment. This makes 
the state of ensembles of cold and/or ultra cold atoms, 
mico/nanomechanical oscillators and the vibrational de- 
gree of freedom of trapped particles, as well as the inter- 
face between such systems, perfectly suited to be tested 
via the phase-space approach described here. In fact, our 
formalism holds the potential to embody a very powerful 
instrument for the analysis of the sharing of quantum 
correlations among the elements of a multipartite hy- 
brid system comprising effective bosonic modes of various 
natures, as well as truly many-body arrays of identical 
trapped particles whose vibrational degrees of freedom 
would embody the register of bosons to study. 

In particular, we can explicitly mention linear [3~4j . as 
well as planar [3|| and multipole ion traps [36j], which 
are able to accommodate ions trapped in unidimensional, 
as well as bidimcnsional configurations (so-called ionic 
crystals: bidimcnsional structures are formed either via 
crystallization of explicit confinement of multi-ion ar- 
rays). In this case, the vibrational degrees of freedom 
of the trapped ions, which are effectively interacting due 
to Coulomb repulsion, give rise to interesting long-range 
couplings that are likely to result in multipartite states 
whose nonlocal character can well be studied by the 
means of our tools (2(| [37| ■ Similar considerations hold 
for ultra cold atoms loading an intracavity double-well 
potential, where in the so-called two-mode approxima- 
tion [38[, a tripartite bosonic state of a cavity field and 
two atomic modes can be established. 

Finally, it will be interesting to address the case of mul- 
tipartite optomechanical devices comprising a mechani- 
cal system interacting with multiple optical modes, and 
hybrid optomechanical systems comprising vibrating me- 
chanical modes, such as the situations addressed explic- 
itly in Refs. |4^[4l|, as well as hybrid configurations in- 
cluding mechanical systems that are coupled to ultracold 
atoms (such as a Bose- Einstein condensate) [331 l42l EU 
or an optically trapped dielectric sphere j44j|. In Fig. [H] 
we show the sketch of some of such instances of physical 
systems, which arc discussed in details in the references 
mentioned in the figure caption, to which we refer. 
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(a) (b) 



(c) 




FIG. 6: Exemplary physical systems where our phase-space 
formalism could be applied to investigate genuine multipartite 
nonlocality: (a) An ultra cold atomic system loaded into an 
intra-cavity double- well potential [3S| ; (b) a Coulomb crystal 
of ions confined on a planar on-chip trap (the dashed tri- 
angle identified a subsystem of three ions whose state can 
be studied with our proposed tools. Red [Blue] arrows show 
nearest-neighbor [next-to-nearest-neighbor] interactions) [3^ |; 
(c) a mechanical membrane in the middle of a cavity [4CJ]; (d) 
a doped micro-toroid coupled to a fibre (4j|; (e) a hybrid 
optomechanical device including a mechanical mode and a 
trapped particle [44|] or Bose-Einstein condensate [43 |. 



VI. CONCLUSIONS 

We have formulated Svetlichny-type Bell tests for s- 
parameterized quasiprobability functions in phase space 
and tested them using interesting three-mode CV states. 
GHZ- and W-type nonlocality can be distinguished, with 
our tools, by properly adjusting the value of s. GHZ-typc 
nonlocality with many photons appears to be witnessed 
using local parity measurements, while single-photon W- 
type nonlocality requires on-off measurements (i.e. the 
Husimi function). This reflects the fact that testing gen- 
uine multipartite nonlocality would require the pondered 
choice of local measurement settings, depending on the 
type of state being studied, in contrast to the violation 
of MK, which occurs for a wide ranges of values of s. 

We observe that not all W-class entanglement can vio- 
late a Svetlichny inequality, indicating that the nonlocal 
character of genuinely multi-mode entangled states docs 
not always coincide with genuine multi-mode nonlocal- 
ity. In studying the behaviour of Svetlichny parameters 
for tripartite entangled states, we found that ECSs are 
useful to demonstrate genuine three-mode nonlocality in 
the whole range of amplitudes of the state components 
and that, the stronger GHZ-typc nonlocality, the more 
fragile it is under detrimental effects. Differently, a pro- 
nounced W-typc nonlocal character is robust to detection 
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inefficiency and local damping. 

Our results consistently show that multi-photon states 
are endowed with very fragile genuinely multipartite non- 
locality content. As our analysis reveals, W-type nonlo- 
cality is exhibited by states of a small amplitude (asymp- 
totically, multimode single-photon states), while GHZ- 
type nonlocality is inherent in GHZ-class entangled states 
of many photons. This explains why the GHZ-type non- 
locality is found to be more prone to environmental ac- 
tions. 

Our phase-space inequalities embody useful and pow- 
erful tools for the investigation of multipartite quantum- 
ness in quantum technology settings of current experi- 
mental interests. As discussed in this paper, there is large 
fan of physical configurations, involving effective bosons 
of various physical embodiments, that can benefit from 
the application of phase-space methods for the tests of 
multipartite nonlocality, a task that we will pursue in 
the close future. 
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Appendix 1 



In this Appendix we report the explicit form of the 
s-paramctcrized MK parameters for the two ranges of 
values of s stated in the main body of the paper. Follow- 
ing the approach described in Sec. HH we obtain the MK 
parameters 



n 3 (l-s) 6 7r 2 (l-s) 4 s 
■M{-i< s <o} = V W3 (at,f3,j;s) H [D W2 (a,/3; s) + % 2 (J3,j; s) + XV 2 (7, a; s) 

+ 7r(l-s) 2 s 2 Yj W 1 (6;s) + 2s 3 , 

7W {s <_ 1} = ^ 3 (1 - sfV W3 (a, (3, 7; s) + tt 2 (1 - s) 2 [XV, (a, (3; s) + XV 2 (/3, 7; s) + XV 3 (7, a; s)] 
+ 2tt(1- S ) W 1 (6;s)-2, 



(A-l) 



where we have introduced the function 



and its marginals 



V W3 (a, f3, 7; s) = W 3 (a, 0, 7'; «) + W 3 (a, /?', 7; s) 
+ W 3 (a',p,r,s)-W 3 (a',p' n ';s) 

(A-2) 

and, say 



V W2 (a, p- s) = W 2 (a, (3; s) + W 2 (a, s) 
+ W 2 (a',j3;s)-W 2 (a',l3';s). 



(A-3) 



which arc defined in terms of the set of measurement 
outcomes 5 = {5, 6'} with S = a,/3,7, the three-mode 
correlator 



W 2 {a,p;s) = J W 3 {a,P, 1] s)d 2 1 , 
Wi(a;s)= ( V 3 (a,/3,7;s)d 2 /3d 2 7. 



Finally, we provide the form of the Svetlichny parame- 
ter expressed in terms of the Wigner and Husimi Q func- 
tions, which read 



W 3 (a,j3, r ,s) 



7T 3 (1 - S) 3 



Tr[pfl(a; s)®IT(/3; s)®n( 7 ; s)] 



S. 



s=0 



- [D W3 (a, (3, 7; 0) ± XV 3 (a, /3, 7; 0)] , (A-4) 
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S.=_i = 8^ 3 [P Q3 (a,/3, 7 ) +% 3 (a,/3, 7 )] 
+ 87r 2 [Q(a,/3') + Q(/3,7') + Q(7,a')] 
+ 87r 2 [Q(a',/3) + 0(/3', 7 ) + Q( 7 ',a)] 
+ 4tt £ [Q(*)+Q(*')], 



(A-5) 



S— a,/3,7 

respectively, where 
O^ 3 (a,/3 )7 ;0) 



(A-6) 



^(a',/3', 7 ) + ^(«',/?,7') 
W(a JJ 8' J 7 , )-W(a JJ 8,7). 

Appendix 2 



In this Appendix we provide the formal details behind 
the replacements to be operated on the s-parameterized 
quasi probability distributions under the effects of ineffi- 
cient detectors and local amplitude damping at non-zero 
temperature. 

We start considering the effects of imperfect detec- 
tion, which can be in general described using a vir- 
tual beam splitter placed in front of a detector with 
perfect efficiency. Detection efficiency is thus associ- 
ated with the transmittivity 77 of the virtual beam split- 
ter, which changes the true photon-number distribution 
P(n) of a signal into the measured quantity P r) (m) = 
Er=m P W(m)(l 7 v) n - m V m @- Correspondingly, 
the measured quasiprobability function at the origin of 
phase-space (i.e. for a — 0) reads 



W v (0;s) 



7r(l - S) ^ VS - 1 



Pr,(m) 



WO; 



I— £—77 
V 



(B-l) 



V 



The s-parameterized quasiprobability function measured 
by a detector with efficiency r\ can thus be identified with 
a rescaled quasiprobability function characterised by the 
inefficiency-dependent parameter [30j 



1 — S — T] 



(B-2) 



Let us now describe the state evolution in local ther- 
mal environments. The scope is to account for such dy- 
namics as the changes of s. The effect of the local ther- 
mal baths can be modelled by the mixture, at a beam 
splitter, of the state of the mode under scrutiny with a 
thermal field. This help modelling the dynamics of the 
three modes in terms of a Fokker-Planck equation for 
the s-parameterized quasiprobability distribution func- 
tion reading [46j |. 



dW(a,p, r ,s;T) _ r 
Ot 2 



as 6 ' 







85* 



+ r, 2 



H 9585 

o— a,p ,7 



E 

(5— a,/3,7 

d 2 

W(a,l3,r,s;T), 



W(a,P,r,s;r) 



(B-3) 



where r is the system-bath interaction time and T is the 
energy decay rate into the environmental baths each of 
average thermal photon number h. The evolution can 
be cast into the form of a convolution between the quasi 
probability function of the three system modes and those 
of the three environmental baths at thermal equilibrium. 
That is 



t ( '( 



^ / d 2 a 'd 2 /3'd 2 7 'n 5 '^ th (<5';s) 
( T ) J 



y[ -a-r(T)^^-r(r)^ i7 -r(r)y 



t(r) 



t(r) 



t(r) 



(B-4) 

where r(r) and t(r) have been introduced in the body of 
the paper, and 



W th (a;s) 



n(2n + 1 - s) 



- e 2A+1- 



(B-5) 



is the s-parameterized quasiprobability function for the 
thermal state of average thermal photon number n. The 
environmental effects can thus be identified with tempo- 
ral changes of the quasiprobability function as stated in 
Eq. (fT8|) . which concludes our formal assessment. 
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